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Abstract 



We investigate the relationship between the Lyapunov exponents of periodic 
trajectories, the average and fluctuations of Lyapunov exponents of ergodic 



(N 
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f^ ■ trajectories, and the ergodic autocorrelation time for the two-dimensional hy 

^. 

^D ■ perbola billiard. We then study the fluctuation properties of the ergodic Lya 

G^ 

- punov spectrum of classical SU(2) gauge theory on a lattice. Our results are 

consistent with the notion that this system is globally hyperbolic. Among the 
D . many powerful theorems applicable to such systems, we discuss one relating 



to the fluctuations in the entropy growth rate. 
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I. INTRODUCTION 

Extensive experimental efforts are under way at Brookhaven National Laboratory and 
CERN to produce and investigate the new deconfined, chirally symmetric high-temperature 
phase of QCD, usually called the quark-gluon-plasma (QGP). While the very high energy 
densities generated in high-energy nuclear collisions virtually guarantee that some new state 
of matter is reached, there are still important unresolved theoretical problems relating to 
the description of this state. One missing, critical ingredient is a non-perturbative approach 
to dynamical QCD processes far from thermodynamical equilibrium. 

The study of non-equilibrium dynamics of relativistic quantum fields is currently an 
active area of theoretical research IQ. Approaches that go beyond perturbation theory 
include descriptions in terms of probabilistic transport equations 0, and deterministic or 
stochastic classical equations for the infrared degrees of freedom of the quantum fields p,H . 
In the special, but important case of non-Abelian gauge theories, the extreme infrared limit 
has been long known to correspond to a djTiamical system exhibiting classical as well as 
quantum chaos f^J^. Several years ago, this result was extended to spatially varying, lattice 
regulated Yang-Mills fields by numerical calculation of the maximal Lyapunov exponents 
and the complete ergodic Lyapunov spectrum of classical SU(2) gauge theory |0-§. The 
most intriguing results with implications for relativistic heavy ion physics are: 

1. The ergodic Lyapunov spectrum looks exactly as expected for a globally hyperbolic 
system. 

2. The largest Lyapunov exponent appears to be related to the plasmon damping rate as 
predicted by high temperature perturbation theory |10 . 



3. The magnitude of the maximal Lyapunov exponent for SU(3) indicates a rapid ther- 
malization of gluons in heavy-ion collisions. 

These results suggest the extension of this approach to a systematic semi-classical description 
of the dynamics of Yang-Mills field theories. The success of such an approach will ultimately 



depend on one's ability to find practical methods for the application of periodic orbit theory 
to systems with many degrees of freedom. We discuss a very first step in this direction. 

We present results of an investigation of the relation between the Lyapunov exponents of 
periodic and ergodic orbits. Periodic orbit theory, in the framework of the thermodynamic 
formalism, makes detailed predictions for the statistical properties of Lyapunov exponents of 
generic orbits in Anosov systems, but few studies of these relations appear to have been made 
for specific chaotic, non-linear dynamical systems. This motivated our numerical study of the 
relation between the Lyapunov exponents of periodic orbits and generic trajectories in a sys- 
tem for which the Lyapunov exponents of periodic orbits (henceforth simply called "periodic 
Lyapunov exponents" ) are known for all orbits below a certain period: the two-dimensional 
hyperbola billiard [jll|. For this system, therefore, powerful mean value theorems can be 



invoked to predict analytical relations which can be checked numerically. Below, we present 
our numerical results confirming the general connection between the Lyapunov exponents 
for ergodic and periodic orbits as well as for their fluctuations. 

We then dicuss the corresponding properties of the Lyapunov exponents of ergodic tra- 
jectories for classical SU(2) Yang-Mills theory on a three-dimensional lattice. The observed 
similarities suggest that this system is also globally hyperbolic and could, in principle, be 
treated within the framework of periodic orbit theory. Our conjecture yields a prediction for 
the fluctuation properties of the ergodic Lyapunov exponents which is verified numerically. 
On the basis of the relation between these fluctuations and the fluctuations of the entropy 
growth rate we obtain a prediction of the magnitude of entropy fluctuations as a function 
of space-time volume. We find that for the conditions occurring in high energy nuclear 
collisions these fluctuations are expected to be very small, in agreement with observations. 

We emphasize that it is presently impossible to predict how far this approach will carry 
toward a description of the dynamics of non-equilibrium processes in QCD. Classical Yang- 
Mills equations can only be used to estimate a very limited number of dynamical parameters 
of the QGP, namely those which have a well-defined classical limit, such as the logarithmic 
entropy growth rate, d(\nS)/dt, but not quantities such as the energy or entropy density. 



Our analysis is only relevant to the fluctuation properties of such, essentially "classical" 
quantities. However, we hasten to point out that, independently of the speciflc application 
considered here, an improved understanding of the connection between quantum field theory 
and periodic orbit theory is of fundamental theoretical relevance for non-linear dynamics in 
general. To our knowledge for the first time, we propose a general relationship between the 
mean periodic Lyapunov exponents of a dynamical system, its mean ergodic Lyapunov ex- 
ponents, and the ergodic autocorrelation time. This general relationship makes it possible to 
extract important new information for any higher- dimensional system for which the explicit 
construction of the periodic orbits is practically not feasible. 

The basic assumption underlying periodic orbit theory is that the periodic orbits sample 
the phase space of a non-linear dynamical system in such a manner that its averaged proper- 
ties can be systematically reconstructed from the properties of the periodic orbits. For each 
such orbit there is a spectrum of characteristic Lyapunov exponents that describe how fast 
the separation between neighboring orbits increases with time. While periodic orbit theory 
is an extremely powerful tool, its range of applicability is strongly limited by the difficulties 
encountered in determining the complete set of periodic orbits. For any field theory with 
its potentially infinitely many degrees of freedom, the task of numerically constructing the 
periodic orbits looks hopeless. It is, however, relatively easy to obtain ergodic Lyapunov 
exponents by numerical integration of the equations of motion 0]. Since it seems plausible 
that every ergodic trajectory eventually comes close to any periodic orbit, any infinite er- 
godic orbit should sample all periodic orbits. Thus it appears as a natural conjecture that 
the average properties of ergodic Lyapunov exponents and the average properties of periodic 
Lyapunov exponents should be related. It is this relationship that we want to discuss in the 
following. 



II. GENERAL RELATIONS 

Before we investigate and confirm the relationship between ergodic and periodic orbits 
for a simple but non-trivial system for which all periodic Lyapunov exponents (up to a 



certain period) are known, namely, the two-dimensional hyperbola billiard [11|, we review 
some general relations between Lyapunov exponents of periodic and generic trajectories. In 
the next section, we will compare these analytic predictions for the properties of the ergodic 
Lyapunov exponents Ar with those obtained by numerical integration of a randomly chosen 
ergodic trajectory x(t) = Xo(t) + 6x(t): 

A,= hm hm-ln^^M, (1) 



where the index r indicates the random starting point. (We remind the reader that for a 
fully ergodic system this yields the maximal ergodic Lyapunov exponent, which for d = 2 
degrees of freedom is the unique positive exponent.) 

In a Hamiltonian hyperbolic dynamical system with d degrees of freedom ergodicity 
implies that the sum of its (i — 1 positive ergodic Lyapunov exponents can also be obtained 
as the ergodic mean of the local expansion rate, 

lim Kit) = lim - /* x{x{t')) dt' = J2 Kj = Jiks ■ (2) 

t—yoo t^oo t Jo ~^ 

Here /iks denotes the Kolomogorov- Sinai entropy and 

\ y ' / expanding 

is the local rate of expansion along the trajectory x{t) . Due to the equidistribution of periodic 
orbits in phase space it is possible to evaluate the ergodic mean in (H) by weighted sums 
over periodic orbits. In fact, for hyperbolic systems the thermodynamic formalism allows to 
express certain invariant measures on phase space in terms of averages over periodic orbits. 



see, e.g., [|T2|,IT3|. One is in particular able to obtain a relation that establishes a direct 



connection between the positive ergodic Lyapunov exponents A^,- and those of periodic 



orbits. Labelling periodic orbits by u, and denoting their periods and positive Lyapunov 
exponents by Ty and X^^j, respectively, this relation reads 



> Ar ,• = hm — ^^ ^-^ ^ -^- — ^ ^ , (4) 



P(/?)=hm^ln ^ exp -/5^A,,,rJ , (5) 



i— >oo 

where e > is arbitrary. Within the thermodynamic formalism the topological pressure 
P{I3) was introduced as a useful tool to analyze invariant measures on phase space in terms 

of periodic orbits as, e.g., in (§). This function can be expressed as 

/ d-i 
In y^ exn I —Fi 

and it is not difficult to derive from (^ that P{I3) is monotonically decreasing and convex. 
The exponential proliferation of the number of periodic orbits immediately implies that 
P(0) = htop (topological entropy). Moreover, the arithmetic average of the sum of the 
positive periodic Lyapunov exponents is given by A = — P'(0). The relation @) then follows 
from (!]) and from the non-trivial identity — P'(l) = /iks- One also concludes that the 
three quantities measuring a mean separation of neighboring trajectories are ordered in the 
following way: A > /i^p ^ ^ks- For further information see, e.g., |]13 . 



Our next goal is to investigate the fluctuations of the local rate of expansion (^, when 
integrated up to a sampling time tg, about its ergodic mean (0). We recall that this quantity 
was denoted as h^{ts) in (^. For (uniformly) hyperbolic dynamical systems one expects 
that observables sampled along ergodic trajectories up to time tg show Gaussian fluctuations 
about their ergodic mean. Indeed, in many cases a central limit theorem holds true that 
also predicts the widths of these Gaussian to scale as t'^^"^ for large sampling times t^. More 
precisely, Waddington [Q has shown that for Anosov systems (i.e., fully hyperbolic systems 
on compact phase spaces) the difference 



%[K{U)-hKs] (6) 

shows Gaussian fluctuations with variance P"(l) in the limit ts -^ oo. This means that 



A/ir(ts) ~ V^"(l)As , ts ^ OO . (7) 



According to (^) the quantity -P"(l) can be expressed in terms of periodic orbit sums as 



P"(l) = lim t 



Eu{EjK,j) exp (- Ej A^jT^) fEu{EjK,j)exp(^-J2jK,jT^ 



Ei.exp(-E,A,.jT, 



Eu exp (- Ej K,jT^ 



(8) 



On the other hand, the variance of the distribution of the periodic Lyapunov exponents is 
related to -P"(0), since 



P"(0) = Um t 



E 



t<T„<t+e 



{J2j=l K,j 



J2t<Tv<t+e (5_/j=l ^u,j 



E 



t<Ty<t+e 
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E 



t<Ty<t+e 



1 



(9) 



For the hyperbola billiard this variance was calculated numerically by Sieber |^, who 
found Gaussian distributions of the positive Lyapunov exponents of periodic orbits with N 
bounces off the boundary. For large A^ the widths of these Gaussian scale like 

0.199 



CTat 



A^ 



(10) 



Taking into account that the mean length of periodic orbits with A^ bounces scales as 



tjv ~ 2.027A^ |jTT|, this yields a prediction for the width of the distribution of periodic 
Lyapunov exponents expressed as a function of t that scales as 

0.283 



AA,(t) 



x/t 



(11) 



in the limit of long periodic orbits. One hence concludes that P"(0) = 0.08. 

The variance of the fluctuations (10) can also be related to the autocorrelation function 



a(r) = (x(x(r))x(f(0)))-(/iKs)' 



(12) 



of the local ergodic Lyapunov exponents, where (...) denotes a phase space average. In 
order to derive this connection one averages the square of (^) over phase space, which then 
leads to 



t{Ahr{t)f 



1 /•+* 



t J-t 



(t- |r|)a(r) dr 



(13) 



A connection with the topological pressure can be established because (|) and (0) imply that 
the autocorrelation function (|T2]) vanishes faster than 1/r as r ^ oo. One can therefore 
perform the limit t ^ oo on both sides of (pTSD, yielding 



/ + 00 
air) dr . (14) 

-oo 

Finally we want to discuss the probability for deviations of the sum of the positive ergodic 
Lyapunov exponents, sampled over time t, from its ergodic mean /iks- To this end let Pt{h) 
denote the probability density for hj.{t) to have a value h. Waddington has shown [|I4| that 
for Anosov systems which are such that P"{P) 7^ for all P, this probability density has the 
form 

pt{h) = f{h) Vi expi-gih)t) , (15) 

where f{h) is a complicated, though uniquely fixed function. Moreover, 

g{h) = mi{hp + P{(3 + l)} (16) 

P 

is a strictly convex, non-negative function with a unique minimum at the ergodic mean 
hmin = ^KS, where 5'(/iks) = 0. This means that for large t the probability of large deviations 
of hy.{t) from the ergodic mean is exponentially small. 

III. THE TWO-DIMENSIONAL HYPERBOLA BILLIARD 

We test the above statements in the two-dimensional hyperbola billiard, for which all 
periodic orbits and their Lyapunov exponents are known up to certain orbit period [|Tl|. In 
order to be able to compare our numerical results with the analytical predictions, which 
are based on periodic orbits in a restricted length range, we have limited the motion into 
the arms of the hyperbola billiard, using the cut-off |a;|, \y\ < xnm = lO/y^ and reflecting 
the motion horizontally or vertically at the boundary. We have not studied the dependence 
of our results on xum in any systematic fashion, but a cursory exploration did not reveal a 
significant dependence. 

8 



Our numerical result for the KS-entropy was obtained as /iks = K = 0.575 by exploiting 
the relation (|1]) for the positive ergodic Lyapunov exponent. In |Tl| the arithmetic average 



of the periodic Lyapunov exponents and the topological entropy have been determined nu- 
merically as A = 0.703 and htop = 0.5925, respectively, so that the ordering A > htop > h^s 
is respected. This provides a non-trivial test since the general theoretical statement has only 
been proven for uniformly hyperbolic systems (Anosov systems) on compact phase spaces. 
The hyperbola billiard is only non-uniformly hyperbolic and, moreover, without the imposed 
cut-off its phase space fails to be compact. 

For the (cut-off) hyperbola billiard we found that the distributions of the ergodic Lya- 
punov exponents (|lD that we determined numerically up to sampling times t^ are very well 
described by Gaussians, see Fig. |], if the sampling time is not too small (tg ^ 1). For 
small sampling times, most of the phase space divergence occurs during intervals t^ when 
the trajectory reflects off the hyperbolic boundary, making the distribution of Kits) = Kits) 
strongly non-Gaussian in the limit t^ —>■ 0. We made power-law fits of the form at^^ to the 
dependence of the widths of these Gaussians on tg. This gave the result (see Fig. 0): 

Ah,{Q^OMt;^/^ . (17) 

We also determined the correlation function a{T) for the hyperbola billiard by sampling 
h^(ts) for in small intervals ts = l/(2v^). The result is shown in Fig. ^. Clearly, a{T) falls 
off rapidly with a time constant of about tc = 6. Therefore, we can test the relation (|13D by 
integrating the right-hand side numerically. For t = 28.3, corresponding to the lower plot 
in Fig. m, we obtain in this way the prediction Ah,. = 0.197 with an estimated numerical 
uncertainty of about 25%. The value obtained from the Gaussian fit to the histogram in 
Fig. |1] is Ahj. = 0.159. The quality of this agreement must be judged with the fact in mind 
that the correlation function a{T), as well as the distribution x{x{t)) are highly singular for 
the hyperbola billiard in the limit r -^ 0. 



IV. THE SU(2) GAUGE THEORY ON A LATTICE 

Let us now turn to a comparison with results obtained for ergodic orbits in the classical 
SU(2) Yang-Mills theory regularized on a lattice. In ^ the complete (positive) Lyapunov 
spectra were obtained for lattice volumes Lr" with L = 1, 2, 3. We have extended these 
calculations to the lattices of size L = 4, 6. All our calculations were performed for an 
average energy per plaquette Ep ^ 1.8. For sufficiently long trajectories and fixed energy 
per lattice site the Lyapunov spectrum has a unique shape, independent of the lattize size, 
as shown in Fig.|^. Indeed, for a completely hyperbolic system, physical intuition requires 
that the Kolmogorov-Sinai entropy — P'(l) is an extensive quantity. For this to be true, 
the sum over all positive Lyapunov exponents must scale like the lattice volume L^ and 
the shape of the distribution of Lyapunov exponents must be independent of L. Figure |^ 
confirms this expectation. 

In Fig. ^ we show distributions of the sum over positive Lyapunov exponents as a function 
of the length of the sampled ergodic trajectories (obtained as function of the sampling time tg 
on a single, very long trajectory). Obviously, the distributions are nicely fitted by Gaussians 
whose widths decrease like l/y^ (see Fig. ^. This behavior is identical to that of the 
two-dimensional hyperbolic system studied before (cf. Fig. ^. We also determined again 
the autocorrelation function a^r) defined in (|T2|) by sampling the distribution pt{h) with 
small time steps (see top part of Fig. ^). For the L = 4 lattice the result is shown in the 
lower part of Fig. ^. This allows us to test the relation ([T3|) connecting the a{T) with the 
variance of the ergodic Lyapunov exponents. Using (|13D we obtain the value Ah^ = 0.88 for 
ts = 6, whereas the Gaussian fit to the sampled distribution shown in the top part of Fig. |^ 
is Ah, = 0.83. 

One can also read off from the distributions shown in Fig. ^ how the widths of the 
Gaussians scale with the lattice size L. To a very good approximation we find that it is 
proportional to -\/Z^. If one includes the sampling time dependence, the variance of /iks 



scales like JL^/ts. As the mean value hxs of the distribution pt{h) scales like L^, this result 
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confirms the Gaussian nature of the fluctuations. Our result also has important consequences 
for heavy-ion collisions. If fluctuations are Gaussian with a dimensional scale given by the 
mean maximal ergodic Lyapunov exponent, which is found numerically to be of order (0.5 
fm)~^ 1^, then for typical volumes and reaction times encountered in nuclear reactions the 



relative fluctuations must be very small, of order w(0.5fm)'^/(5fm)'^ = 0.01. This result is in 
agreement with a recent measurement of the fluctuations in relativistic heavy-ion collisions, 
which show that the primary event-by-event fluctuations in the mean value of the transverse 
momentum do not exceed 1 percent [^ . 



Let us stress that while it is consistent to assume that the SU(2) gauge theory treated as 
a classical field theory on the lattice is a hyperbolic system, our positive evidence is limited. 
It should be clear that it is impossible to exclude, by numerical calculations for a limited 
number of trajectories, that there are regions in the high-dimensional phase space of our 
lattice field theory which are not hyperbolic. (Then the SU(2) field on the lattice would not 
be an Anosov system.) Also it is unproven, though highly probable, that the addition of the 
quarks will not change the picture. 

V. CONCLUSIONS 

We have shown by numerical simulations that for a two-dimensional billiard the mean 
values for the ergodic and periodic Lyapunov exponents and their fiuctuations as a function of 
trajectory length (i.e. time) are closely related. We have derived a general relation between 
their mean values and checked it numerically. This demonstrates that we understand the 
relationship between ergodic and periodic Lyapunov exponents for the hyperbola billiard. 
We have than analyzed in a similar way classical SU(2) gauge theory on a lattice. For 
all investigated properties we found good agreement with the expectations for a globally 
hyperbolic (Anosov) system. We conclude that for all quantities of interest which have a 
well-defined classical limit (like the growth rate of entropy after the initial energy deposition 
by hard interactions) the probability for large fiuctuations should be exponentially small. 
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For typical high-energy heavy-ion coUisions (Pb-|-Pb) such fluctuations are estimated to be 
at most of the order of a few percent. 
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FIG. 1. Distribution the calculated local Lyapunov exponents for ergodic trajectories of two 
different length ts in the hyperbola billiard. 
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FIG. 2. The widths of the Gaussians ihustrated in Fig. IT] as a function of tg together with a fits 
of the form atg (solid hne) and at~^ with b = 0.45 (dashed hne). 
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FIG. 3. Temporal autocorrelation for the local Lyapunov exponents determined along an er- 
godic trajectory in the two-dimensional hyperbola billiard. The dashed line is a fit of the form 
aexp(— r^/tg) yielding tc ~ 6. 
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FIG. 4. Distribution of numerically obtained ergodic Lyapunov exponents for a classical SU(2) 
gauge theory on lattices of size L = 2,4, 6. The index i numbers the Lyapunov exponents and the 
abscissa is scaled with L^. 
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FIG. 5. Distributions of the sum of the positive ergodic Lyapunov exponents for different values 
of the trajectory length ts for classical lattice SU(2) gauge theory on lattices of size L = 4, 6. 
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FIG. 6. The widths of the Gaussians illustrated in Fig. |^ and scaled with {ts/L^y'^, as a 
function of tg- 
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FIG. 7. Top: The distribution of the sum of local expansion rates hj- for L=4 and a short 
sampling time tg = 0.1, together with a Gaussian fit. Because the globally expanding phase space 
volume may be locally contracting, the distribution has a small tail extending to negative values 
of h^. This tail disappears (becomes exponentially small) for large values of tg- Bottom: The 
autocorrelation function a[T) for this distribution. 
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